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We study the consequences of supersymmetry for primordial black hole (PBH) abundance con¬ 
straints. PBHs will emit supersymmetric particles throughout their evaporation if their mass is less 
than about lO^^g. In most models of supersymmetry the lightest of these particles, the lightest 
supersymmetric particle (LSP), is stable and will hence survive to the present day. We calculate 
the limit on the initial abundance of PBHs from the requirement that the present day LSP density 
is less than the critical density. We apply this limit, along with those previously obtained from 
the effects of PBH evaporation on nucleosynthesis and the present day density of PBHs, to PBHs 
formed from the collapse of inflationary density perturbations in the context of supersymmetric 
inflation models. If the reheat temperature after inflation is low, so as to avoid the overproduction 
of gravitinos and moduli then the lightest PBHs which are produced in significant numbers will be 
evaporating around the present day and there are therefore no constraints from the effects of the 
evaporation products on nucleosynthesis or from the production of LSPs. We then examine mod¬ 
els with a high reheat temperature and a subsequent period of thermal inflation. In these models 
avoiding the overproduction of LSPs limits the abundance of low mass PBHs which were previously 
unconstrained. Throughout we incorporate the production, at fixed time, of PBHs with a range of 
masses, which occurs when critical collapse is taken into account. 
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I. INTRODUCTION 

Primordial black holes (PBHs) may form in the early 
universe via a number of mechanisms the simplest 
of which is the collapse of large density perturbations [Q . 
Due to quantum effects PBHs evaporate, mimicking the 
emission from a black body with finite size and temper¬ 
ature Tbh where 

and Mbh is the PBH mass in grams. The standard pic¬ 
ture of PBH evaporation is that all particles which appear 
elementary at the energy scale of the PBH and have rest 
mass less than Tbh are emitted directly ^ . For instance 
PBHs with Tbh above the QCD quark-hadron transition 
scale Aqh ~ 250 — 300 MeV, emit relativistic quark and 
gluon jets which then fragment into photons, leptons and 
hadrons as in high energy accelerator collisions. 

The mass loss, in grams per second, is given by 

= -5.34 X 10"V(Mbh)Mb^ , (2) 

where /(Mbh) is a function of the number of species 
emitted and is normalised to 1 for PBHs with Mbh ^ 
lO^^gwhich emit only massless particles (see ref. and 
Sec. H for more details). The effects of the evaporation 
products allow constraints to be put on the initial abun¬ 
dance of PBHs over a range of masses. PBHs with mass 
in the range 10®g < Mbh < lO^^g, would have evapo¬ 
rated after nucleosynthesis and could have a number of 


effects on the successful predictions of nucleosynthesis Q| 
whilst PBHs with Mbh ~ 5 x 10^‘^g are evaporating to¬ 
day and the number density of photons produced must 
not exceed the observed y-ray background . In the case 
of lighter PBHs with Tbh > 100 GeV the fundamental 
particles emitted, and hence /(A/bh)) depend upon the 
particle physics model assumed. 

Gurrently the widely accepted extension of the Stan¬ 
dard Model is supersymmetry [^, where each Standard 
Model particle has a supersymmetric partner known as a 
sparticle. Motivated mainly as an attempt to understand 
why the weak scale is much smaller than the Planck scale 
(known as the gauge hierarchy problem 0), supersym¬ 
metry also leads to the unification of gauge couplings at 
an energy of about 10^® GeV [||. The phenomenology 
of the sparticles is complicated being governed by up to 
105 independent and unknown parameters. In the sim¬ 
plest models the sparticles have masses of order 100 GeV 
and there is a multiplicatively-conserved quantum num¬ 
ber known as R-parity, where Standard Model particles 
have R= -1-1 and sparticles have R= —1. Gonsequentially 
heavier sparticles decay into lighter sparticles and the 
Lightest Supersymmetric Particle (LSP) is stable, since 
it has no allowed decay mode. In most models the LSP is 
non-relativistic at freeze out and is therefore a candidate 
for the cold dark matter (GDM) 


*In gauge-mediated models of supersymmetry breaking the 
LSP could be a gravitino with mass ~ 1 keV which would 
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The LSP must be neutral and weakly interacting 0 
as otherwise it would have condensed, along with the 
baryonic matter, into astrophysical structures and the 
resultant abundance of anomalous heavy isotopes would 
exceed observational limits ||^ . The LSP may therefore 
be a sneutrino or the gravitino but in most supersym¬ 
metric theories it is the lightest neutralino x, which is a 
mix of the supersymmetric partners of the photon, the 
Z boson and the neutral Higgs boson. Throughout this 
paper we therefore use y to denote the LSP. 

Purely experimental searches at LEP have led to the 
limit > 30 GeV [|^. This limit can be tightened 
to m-)^ > 42 GeV by making various theoretical assump¬ 
tions and requiring that the present day LSP density lies 
in the interesting range for GDM: 0.1 < < 0-3 [01 • 

The LSP comoving number density has remained con¬ 
stant since annihilations ceased at the freeze out temper¬ 
ature Tf. This leads to a simple estimate (see Ref. 0 
for a review of this and more detailed calculations) of the 
current LSP density: 




10-3 

(<^ann (xx)i^x)ToWpi ’ 


( 3 ) 


where (CTann(xx)%) is the thermally averaged annihila¬ 
tion cross section. As increases (o'ann(xx)^x) de¬ 
creases so that increases. This leads to an upper 

limit of rrix < 300 GeV [^, or < 600GeV if co¬ 
annihilations with the stau slepton, which are important 
in some regions of parameter space, are taken into ac¬ 
count ]l^ . 

Supersymmetric particles which are produced by the 
evaporation of PBHs after the temperature of the uni¬ 
verse has fallen below Tf will not be able to equilibrate 
due to the inefficiency of annihilations. In Sec. || we cal¬ 
culate the mass range of PBHs which are heavy enough 
to evaporate after LSP freeze out but light enough to 
produce sparticles. The LSPs produced by the decay of 
these sparticles, along with the LSPs evaporated directly, 
will therefore provide an additional contribution to the 
present day density of LSPs. Whilst it is possible that the 
number density of LSPs produced via PBH evaporation 
may be comparable to the freeze out number density, re¬ 
ducing the upper limit on and hence rrix, this would 
require extreme fine tuning of the initial abundance of 
PBHs. We therefore use the conservative requirement 
that the present day density of LSPs produced via PBH 
evaporation after Tf is less than the critical density, to 
constrain the initial abundance of LSP producing PBHs. 

In Sec. HI we outline the resultant constraints on the 
mass fraction of the universe in PBHs formed from the 
collapse of inflationary density perturbations, taking into 
account the formation, due to critical collapse |17|, of 
PBHs with a range of masses at fixed horizon mass. 


Since we are assuming that supersymmetry is the cor¬ 
rect model of particle physics above ~ lOOGeV we must 
apply these constraints in the context of supersymmetric 
inflation models. In Sec. 0 we discuss the constraints on 
successful supersymmetric inflation models and calculate 
and review PBH abundance constraints for two classes 
of inflation model; firstly those where the reheat temper¬ 
ature after inflation is low, and then those with a high 
reheat temperature and a subsequent period of thermal 
inflation. 


II. ABUNDANCE OF LSPS EMITTED BY PBHS 


A PBH will emit LSPs if Tbh is greater than rrix so 
that, using Eq. (0), if Mbh < Ml where 


Ml = 


1.06 X IQi^GeV 


rriy 


( 4 ) 


then the PBH will emit LSPs throughout its evapora¬ 
tion. A PBH with mass greater than Mi will emit LSPs 
during the later stages of its evaporation once its mass 
has fallen below Mi. The the fraction of the PBH energy 
density eventually in LSPs will be reduced by a factor of 
Ml/M bh relative to that for PBHs with Mbh < Mi. In 
this paper we will focus on the case of PBHs formed from 
the collapse of large inflationary density perturbations. 
If the horizon mass at the time the PBHs form is Mh 
then if Mh < Mi the LSP emission will be dominated by 
PBHs with Mbh < Mi. On extrapolating the constraint 
on t he ini tial mass fraction in PBHs which is found in 
Sec. IV B (and incorporating the additional weakening 
factor of Mbh/M l) we can see that for Mh > Mi the 
constraint from LSP emission will be weaker than those 
from the effects of PBH evaporation on nucleosynthesis. 
Therefore we will only calculate the constraints from LSP 
emission by PBHs with Mbh < Mi. 

The PBH lifetime, in grams per second, is ra 


6,24xlO-”Mj„ 

- 70^5 — 


( 5 ) 


so that the temperature at evaporation, Tevap = 
Tpi(tpi/T(MBH))^^^, is given by 


T — 

evap — 




-^BH 


( 6 ) 


For evaporation to occur after LSP freeze out (Tevap < 
Tf), Mbh must be greater than M 2 where, using the fact 
that Tf ~ to^/25, 


M 2 


rriT 


constitute warm dark matter j^. 
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A PBH will therefore evaporate after LSP freeze out and 
also emit LSPs throughout its evaporation if its mass is 
in the range M 2 < Mbh < Mi. 











If at some initial time ti the fraction of the total energy 
density of the universe in PBHs which evaporate after 
freeze out producing LSPs is = Ppbh then 

immediately before the PBHs evaporate 


f Pphh 
\ Prad 

If the PBHs dominate the energy density of the uni¬ 
verse before they evaporate then, since the radiation 
emitted by the PBH dominates the background radia¬ 
tion energy density, soon after evaporation the LSPs will 
come to dominate the energy density of the universe soon 
afterwards. To avoid this we require 


evap 1 ~ /^X.i ^evap 


( 8 ) 


A 

i-A 


< 


evap 


(9) 


The LSPs emitted will initially be relativistic with 
mean energy per particle 3 Tbh The ratio of the energy 
density in LSPs to that in radiation is therefore constant 


Px \ ^ 

Prad / 


'X > 


( 10 ) 


where is the fraction of the PBH mass energy evapo¬ 
rated into LSPs, until the LSPs become non-relativistic 
at Tnr where, since p for relativistic fluids, 


Tnr - Tevap ( 3 ) 


1/4 


( 11 ) 


At any subsequent epoch, with temperature T, before 
matter-radiation equality 


( \ _ ( Px \ 

\ Prad / \ Prad / T 


= [—) 

\ Prad / 


% 


evap 


f ^x \ 


1/4 


evap 


T \3TbhJ 


( 12 ) 


The fraction of the energy density of the universe in LSPs 
at the present day is therefore given by 


Hx,o — llx,®a — ^ 


= 2e. 


^ Tevap i 



^ T„ 1 

\ 3Tbh y 

1 

Px,i Ti 

( rrix ' 


~ /^Xl Teq 

V 3Tbh ^ 

) 


(“) 

\ Prad / , 


(13) 


where ‘eq’ denotes the epoch of matter-radiation equality. 
This relation can be inverted simply to obtain the con¬ 
straint on from the requirement < 1, as a func¬ 
tion of Ti. This constraint is independent of the mecha¬ 
nism of PBH formation and whilst we have assumed that 
all PBHs form at the same time ti it would be simple to 
recalculate the constraint allowing for PBH formation at 
a range of times. 

The relevant mass range and the fraction of the mass 
of the PBH which is evaporated into LSPs both depend 


on /(Mbh)- To fully calculate /(Mbh) as a function of 
Mbh we would need to know the full mass spectrum of 
the sparticles. We can however calculate /(Mbh) in two 
limiting cases. If Tbh then the LSP will be the only 

sparticle emitted (along with photons, gravitons, gluons, 
the three lepton families and six quark flavours) whilst if 
Tbh then all the standard model particles (includ¬ 

ing the Higgs, Z and W bosons) and their supersymmetric 
partners will be emitted. The relativistic contributions 
to /(Mbh) per particle degree of freedom are [^: 

fs^o = 0.267, /,=! = 0.060, 

/,=3/2 = 0.020, /.=2 = 0.007, 

/s=i /2 = 0.147 uncharged, 

fs=i /2 = 0.142 electrically charged. (14) 


A. Tbh ~ rn^ 

In this regime a PBH can emit the 3 lepton families, 6 
quark flavours]^, the graviton, the gluon and the LSP. The 
photon has 2 possible polarisation states and the gluon 
has 48 degrees of freedom giving a total of 50 s = 1 de¬ 
grees of freedom. Each type of neutrino has 2 degrees of 
freedom resulting in 6 neutral s = 1/2 degrees of free¬ 
dom. The electron, muon and tau leptons each have 4 
degrees of freedom, whilst each quark has 12 degrees of 
freedom resulting in a total of 84 charged s = 1/2 states. 
Finally the graviton has s = 2 and 4 degrees of free¬ 
dom. The total contribution of standard model states 
to /(Mbh) is therefore (50 x 0.06) -I- (6 x 0.147) -I- (84 x 
0.142) -I- (4 X 0.007) = 15.84. Finally, as outlined in the 
introduction, the LSP is most likely to be the lightest 
neutralino which is uncharged and has s=l/2, giving a 
final value of/(M bh) = 15.84-^(2x0.147) = 16.13. The 
fraction of the total mass which is evaporated into LSP 
is = 0.294/16.13 = 0.018. 


B. Tbh rn^ 

If Tbh is much larger than all the standard model 
particles and their supersymmetric partners will be emit¬ 
ted. In the Minimal Supersymmetric Standard Model 
there are 2 Higgs doublet fields which give rise to 5 phys¬ 
ical states: the charged Higgs, the light scalar Higgs, 
the heavy scalar Higgs and the pseudoscalar Higgs (see 
e.g. Ref iQ for details). Along with the W and Z bosons 


Mhe LSP mass may be smaller than that of the top quark 
in which case there may be a narrow range of PBH masses 
where the LSP is emitted but the top quark is not, which 
would decrease the value of /(Mbh) calculated here by about 
1.704. 
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these states lead to 9 s = 1 and 5 s = 0 degrees of free¬ 
dom in addition to the s tandard model degrees of freedom 
considered in Sec. II A, so that the total contribution of 
standard model states to /(Mbh) is 17.10 

Each standard model (SM) degree of freedom had a 
supersymmetric (SUSY) degree of freedom with s = 
1/2, 0,1/2 and 3/2 for standard model states with s = 
0,1/2,1 and 2 respectively. This leads to a grand total 
of 95 s = 0 (5 SM, 90 SUSY), 92 charged s = 1/2 (84 
SM, 8 SUSY), 30 uncharged s = 1/2 (6 SM, 56 SUSY), 
59 s = 1 (all SM), 4 s = 3/2 (all SUSY) and 4 s = 2 
(all SM) degrees of freedom, giving /(Mbh) = 46.79. 
The fraction of the total mass evaporated into super- 
symmetric particles is esuSY = 29.69/46.79 = 0.63, with 
= 0.294/46.79 = 0.006 of the total mass evaporated 
directly into LSPs. Whilst the fraction of the PBH mass 
evaporated directly into LSPs is small all the other su¬ 
persymmetric particles emitted will decay rapidly pro¬ 
ducing at least one additional LSPs each so that we can 
set = esuSY when estimate the fraction of the PBH 
mass ending up in the form of LSPs. 


III. PBH FORMATION FROM DENSITY 
PERTURBATIONS 


P{k) = (l^kP) and n is the spectral index, cr^(MH) oc 
during radiation domination ||^. The forma¬ 
tion of PBHs on a range of scales has recently been stud¬ 
ied |p 0 | for power law power spectra and for flat spectra 
with a spike on a given scale. In both cases it was found 
that, in the limit where the number of PBHs formed is 
small enough to satisfy the observational constraints on 
their abundance at evaporation and at the present day, 
it can be assumed that all the PBHs form at a single 
horizon mass. In particular if the power spectrum is a 
power law with n > 1 , as is the case in tree-level hybrid 
inflation models, all PBHs form at the smallest horizon 
scale immediately after re-heating. 

The initial mass fraction of the universe in PBHs with 
masses in the range M 2 < Mbh < Mi, which evaporate 
after LSP freeze out and produce LSPs throughout their 
evaporation, is given by 




/ ' ^p(5(MH))d5(MH) 

Js2 ^ V 


dSiMa), (17) 


where, using eq. 


In order for a PBH to be formed, a collapsing region 
must be large enough to overcome the pressure force re¬ 
sisting its collapse as it falls within its Schwarzschild ra¬ 
dius. This occurs if the perturbation is bigger than a 
critical size Sc at horizon crossing. There also is an up¬ 
per limit of 5 < 1 since a perturbation which exceeded 
this value would correspond initially to a separate closed 
universe | ^. Analytic calculations [Q find Sc ^ 1/3 and 
assume that all PBHs have mass roughly equal to the 
horizon mass at the time they form, independent of the 
size of the perturbation. Recent studies |]^ of the evo¬ 
lution of density perturbations have found that the mass 
of the PBH formed in fact depends on the size of the 
perturbation: 

Mbk = kMu{S — Sc)^ , (15) 

where 7 « 0.37 and k and Sc are constant for a given 
perturbation shape (for Mexican Hat shaped fluctuations 
k = 2.85 and Sc = 0.67). 

In order to determine the number of PBHs formed on a 
given scale we must smooth the density distribution using 
a window function, W{kR). For Gaussian distributed 
fluctuations the probability distribution of the smoothed 
density field p{S{Mh) is given by 

p{SiMB))dS{MB) = 

y/27Ta{Mii) 

where ct(Mh) is the mass variance evaluated at hori¬ 
zon crossing. For power law spectra, P{k) oc fc", where 


Eq.(|^ can be used to translate the constraints on 
into constraints on MMh), which in turn can be used to 
constrain n We can also find the maximum 

fraction of the universe in PBHs of all mass, /Sp. 


A = 


k{S-Scr 

'Sa 


exp - 


2a2(MH) 


dSiMu), (19) 


IV. CONSTRAINTS ON SUPERSYMMETRIC 
INFLATION MODELS 

Since we are assuming that supersymmetry is the cor¬ 
rect model of particle physics above ^ lOOGeV then we 
must apply the constraints on the abundance of PBHs 
formed from the collapse of large inflationary density per¬ 
turbation to models of inflation constructed in the con¬ 
text of supersymmetry. Avoiding the overproduction of 
various relic particles, which would alter the subsequent 
evolution of the universe and wreck the successful pre¬ 
dictions of the standard hot big bang model, leads to 
constraints on supersymmetric inflation models. 

The gravitino (the supersymmetric partner of the 
graviton) has only gravitational interactions and mass of 
order 100 GeV, and will decay after nucleosynthesis. The 
requirement that its decay products do not destroy the 
successful predictions of big-bang nucleosynthesis places 
limits on its abundance. Since the number density of 
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gravitinos is proportional to the reheat temperature af¬ 
ter inflation Trh, this leads to an upper limit on Trh [H) . 
The exact limit depends on the gravitino mass, the avail¬ 
able decay channels and the baryon-to-photon ratio be¬ 
fore the gravitino decays but, very conservatively, 
Trh must be less than 10® GeV. Similarly in almost all 
theories in which supersymmetry is broken at an inter¬ 
mediate scale there are scalar fields, known as moduli 
which typically have the same mass and lifetime as the 
gravitino j^ ]. Avoiding their production after inflation 
requires Trh to be less than 10^® GeV. 

There are several ways to avoid the gravitino and mod¬ 
uli problems. Firstly inflation models can be constructed 
where inflation occurs at a low energy scale so that 
the reheat temperature is automatically low enough to 
avoid these problems. However these models in general 
require fine-tuning . Secondly the reheat temperature 
can be sufficiently low if the inflaton is long lived. For 
instance a model has been constructed [^, using a sin¬ 
glet field in a hidden sector, where inflation occurs at the 
scale of the spontaneous breaking of the gauge symmetry 
which is of order lO^^GeV. In this case, since the infla¬ 
ton has only gravitational strength couplings, the reheat 
temperature is only of order 10® GeV. 

There is however a generic solution to these problems 
in high energy scale inflation models which relies on the 
properties of flaton fields which also arise naturally in su¬ 
persymmetric theories Flaton fields, have vacuum 

expectation values M 10® GeV, even though their 
mass m is only of order the supersymmetry scale, so that 
their potential is almost flat. In the early universe these 
fields are held at zero by finite temperature effects, with 
false vacuum energy density Vq ^ Once the tem¬ 

perature falls below ■, the false vacuum energy den¬ 
sity dominates the thermal energy density of the universe 
and begins to drive a period of inflation known as thermal 
inflation. This inflation continues until the temperature 
drops to T ~ m, at which point thermal effects are no 
longer strong enough to anchor the flaton in the false 
vacuum. Taking M ^ 10®^ GeV gives ^ 10^ GeV 
so that around ln(10^/10®) ^ 10 e-foldings of thermal 
inflation occur, sufficient to dilute the moduli and grav¬ 
itinos existing before thermal inflation but small enough 
to not affect the density perturbations generated during 
the first period of inflation. 

We examine two classes of inflation model; those with 
a low reheat temperature and those with a high reheat 
temperature and a subsequent period of thermal infla¬ 
tion. 


* Specific examples are the dilaton of string theory and the 
massless gauge singlets of string compactifications or, in gen¬ 
eral, any gauge singlet field responsible for SUSY breaking 



FIG. 1. The constraints on the initial mass fraction of 
PBHs, /3i! in inflation models with Trh < 10® GeV . 


A. Low-reheat temperature 


If Trh = 10® GeV, so as to just satisfy the gravitino 
constraint, then, using the relationship between horizon 
mass and temperature in a radiation dominated universe: 


Mh ~ 10 ®® g 


/10®' GeV 

V T 


2 


( 20 ) 


the maximum horizon mass is 1 x 10®^g. Although when 
critical collapse is taken into account PBHs with a range 
of masses are formed at fixed horizon mass, the vast 
majority of PBHs have mass within an order of mag¬ 
nitude of Mr. From eq. must be less than 

5 X 10“®®'. However, for = 30 GeV, only PBHs with 
Mbh < Ml = 3.5 X 10®®g ^ 0.0035A/h are hot enough 
to emit LSPs. The fraction of PBHs which have mass 
this much smaller than the horizon mass is negligible so 
that there is no resultant constraint on f3i. Similarly the 
constraints from the effect of the products of PBH evap¬ 
oration on nucleosynthesis only hold for Mbh < 10 ®®g. 
Therefore in a low reheat temperature inflation model 
only the constraint from the present day density of PBHs 
and, if Trh > 5 x 10® GeV, that on the abundance of 
PBHs evaporating today hold. The constraints on /3i are 
shown in figj^ The ‘U’-shaped dip arises from the limit 
on the abundance of PBHs with mass A/rh 5 x 10®"®g 
which are evaporating at the present day (as calculated 
by Yokoyama and the straight line from the limit 
on the present day density of PBHs which have not evap¬ 
orated. 

The constraints on f3i can be translated into limits on 
the spectral index of the density perturbations |pl| , p^ , p^ 
using Eq. (0) and the scale dependence of a {Mu): 


(t(Mh) 


a{Mo) 


MeJ \Mo ) 


( 21 ) 


where Mq and Meq are the horizon masses at the present 
day and at matter radiation equality and, using the 
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FIG. 2. The constraints on the initial mass fraction of 
PBHs from the present day density of LSPs, in supersym¬ 
metric inflation models with a high reheat temperature and 
a subsequent period of thermal inflation for = 0.018. The 
solid lines shows the constraints on /3x,ii the fraction of the 
universe in PBHs which evaporate after freeze-out and pro¬ 
duce LSPs throughout their evaporation, for (from bottom 
to top) rrix = 600,300, 45 and 30 GeV. The dotted, short 
dashed, long dashed and dot-dashed lines show the maximum 
allowed fraction of the universe in PBHs of all masses, /3i, for 
nix = 600, 300, 45 and 30 GeV respectively. 


COBE normalisation, a{MQ = 10®®g) = 9.5 x 10“®. In 
models with Trh < 10® GeV, the tightest limit on n 
(~ 1.28) arises from the constraint on PBHs which are 
evaporating today. This is weaker than the limits found 
in Ref. firstly because the tightest limits arise from 
the constraints on PBHs of mass Mbh ^ which 

are not formed in low reheat temperature inflation mod¬ 
els. Also the analytically derived value of Sc 1/3) 
used in those calculations is roughly half the new value 
Sc = 0.67, determined from numerical simulations 0. 
which is used in this paper. 


B. High-reheat temperature with thermal inflation 
occurring 

The effects of a period of thermal inflation on the con¬ 
straints on PBH abundance were studied in Ref. . It 
is assumed that thermal inflation commences at T = 10^ 
GeV and continues until T = 10^ GeV, when the fla- 
ton field rolls to its true vacuum state. Assuming for 
simplicity that reheating is efficient then the universe is 
reheated to T = 10^ GeV with the subsequent evolution 
of the universe having its standard form. The duration 
of thermal inflation is negligible compared to the PBH 
lifetime so that its main effect on the PBHs is to dilute 
their density by a factor pi/pf = (of/oi)^ ~ (10^)^ so 
that the constraints on /3i are weakened by a factor of 
^ 1012 . 

If thermal inflation occurs then the right hand sides of 



FIG. 3. A compilation of the constraints on the initial 
mass fraction of PBHs in supersymmetric inflation models 
with a high reheat temperature and a period of thermal infla¬ 
tion. The dotted show the constraints from the present day 
density of LSPs. From left to right the solid lines show the 
constraints from: the effects of PBH evaporation after nucle¬ 
osynthesis, the abundance of PBHs evaporating at present, 
the present day density of PBHs formed before and (on the 
far right) after thermal inflation. 


eqs.(|-|T|) are each multiplied by a factor of 10 i® due to 
the dilution of the PBHs during thermal inflation. Using 
the relation between temperature and horizon mass Tj = 
Tpi(mpi/MH)i^2 the requirement that the PBHs do not 
dominate the universe at evaporation becomes 


< 6.7 X 10® 


/ mh 


( 22 ) 


and the limit from the present day density of LSPs be¬ 
comes 


/3x,i ^ 5 X IQii Teg / 3 X IQi^ y/^ 

l-/3x.i ex Tpi V™Pi/ [mx-MBH/ 

(23) 

To evaluate these constraints on /3x,i we neglect the 
spread in PBHs masses, since the vast majority of PBHs 
will have the same mass to within a factor of a few. If 
Mh is between M 2 and Mi then most of the PBHs will 
have Mbh ~ Afn, whilst if Mh < M 2 the vast majority 
of the PBHs which evaporate after freeze out will have 
Mbh ~ M 2 and similarly if Mh > Mi most of the LSP 
producing PBHs will have Mbh ATi . The tightest con¬ 
straint arises from the present day abundance of LSPs 
for Mh < Ml and also for Mh > Mi if Mi < 2 x lO^^g 
(which is the case for < 30GeV as found experimen¬ 
tally). We calculated the resulting constraints on cr(MH), 
and hence /3i, for 4 sample values of the LSP mass: = 

30, 45, 300 and 600 GeV, using ^ 0.018. As the hori¬ 
zon mass decreases the fraction of the total number of 
PBHs formed which are heavy enough to evaporate after 
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freeze-out decreases. Similarly as the horizon masses in¬ 
creases the fraction of PBHs which are light enough to 
emit LSPs decreases. Emission of LSPs after freeze-out 
constrains j3i for 10®g < Mh < 10^°g. The constraints 
on /3^_i and j3i are shown in Fig. In Fig. || the con¬ 
straints on /3i for = 45GeV are shown, along with 
those from the effects of PBH evaporation on the prod¬ 
ucts of nucleosynthesis ||, the abundance of PBHs evapo¬ 
rating at present and the present day density of PBHs as 
calculated in refs. [|^ and |^. The missing mass range 
corresponds to comoving scales which enter the horizon 
before thermal inflation, and are then pulled back out¬ 
side again during thermal inflation. Any new density 
perturbations are expected to be small, since the energy 
scale of thermal inflation is much lower than the origi¬ 
nal inflationary period, and hence unable to form black 
holes when they re-enter the horizon again after thermal 
inflation. 

If thermal inflation occurs then the temperature at 
which a given comoving scale crosses the Hubble radius 
is changed, whilst the relation between horizon mass and 
temperature remains the same so that the scale depen¬ 
dence of cr(MH) becomes [p^ 


^(Mh) 


a{Mo) 



1 - 

(M 


Meq 


(l-n)/4 


Mo) 


(24) 


The limits on n from the constraints on /3i due to LSP 
emission range from n < 1.36 at Mh ~ 10^^ g to n < 1.33 
at Mh ~ 5 X lO^g. These limits are slightly tighter than 
those from the nucleosynthesis constraints (n ~ 1.34 to 
1.37 if the accurate value of Sc = 0.67 is used). 


V. CONCLUSIONS 

We have examined the consequences of supersymme¬ 
try for PBH abundance constraints. PBHs with mass 
in the range lO^g < Mbh < lO^^g evaporate after LSP 
freeze out and produce LSPs (and other supersymmetric 
particles) throughout their evaporation. In most models 
of supersymmetry the LSP is stable and the requirement 
that the present day density of LSPs does not exceed the 
critical density places a limit on the initial abundance of 
PBHs in this mass range. 

We have studied the constraints on PBH abundance for 
two classes of supersymmetric inflation model; those with 
a low reheat temperature and those with a high reheat 
temperature and a subsequent period of thermal infla¬ 
tion. If the reheat temperature is low the lightest PBHs 


^The sharply rising line from Mh ~ lO^^g arises due to the 
constraint from entropy production after nucleosynthesis on 
PBHs with mass Mbh < lO^^g. 


which can be produced will be evaporating at the present 
day and the constraints from the present day density of 
LSPs and the effects of evaporation on the products of 
nucleosynthesis, which provide the tightest limits on n, 
do not apply. For models with a high reheat temper¬ 
ature and a subsequent period of thermal inflation the 
constraints from the present day density of LSPs pro¬ 
vides the tightest limit on n. 
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